1 

REPORT  DOCUMENTATION  PAGE 

-or-n  JCC';-.  f-J 

0MB  Vo  J'O-i-CISS 

=  c  -ocr-  -c  rur:*?''  •:"  •■  S  v-'-^rtion  :♦  -•;'~ai.O''  ^  *2  '  -^^uf  -“soc'*sc  rc:-jaipg  :r*  :,.-‘.e  -Cf  ••  -  ‘ :  "  *  v  - :  :  ,•  i  ,2  .•  ><  1 

3.“.-  ;  :i'. 3  '>rro«  i'‘a  :  :‘»‘0'e’«f'c  ir-a  “2  :*••  ;iH*'r;icn  '•*  ■}'"*aticn  '>«;no  r*gd--ci!~G  ■•*  s  c--*  2*-^  ■?'. ,  j*®  ;  •  ,  - ,  isc-«  :t  :*  '  >  j 

-t  «•  suq<5«i'C''s  ■  :^ ‘rcucirc  *““%  ■'  /.  Se'^vices.  O'ffcof 3te  *0'’  c.-  ,>?'3.  C“s  ?"•<</%.- 

^3  s 3/  -•’■“g'.cn  ;i  4ria:'i:“*  '•**  •?  :••*.' ire  r-cc*?*  =  ,c>e''-/vC'’<  SeouCG'cr  :  if 't:  1C 'C  J  0 ' 33^  'V^sr -g:  | 

1  AGENCY  USE  ONLY  (Leave  blank)  2. -REPORT  DATE  3._REP^T  TYPE 

December  22,  199  3  j  02'^,*^^!“  ^ 

7  DATES  COVERE/) 

7O  Mol)  93  j  INAL  /' 

4.  TITLE  AND  SUBTITLE 

Nonlinear  Wave  Propagation^^^'^ 

5.  FUNDING  NUMBE^T 

AFOSR  90-0039 

61102F 

2304/BS  i 

R-TR- 9 -  0  0  1 

6.  AUTHOR(S) 

Mark  J.  Ablowitz  AFOw 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

University  of  Colorado 

Program  in  Applied  Mathematics 

Campus  Box  526 

Boulder,  CO  80309-0526 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

1534628 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Air  Force  Office  of  Scientific  Research 

Attn:  Arje  Nachman  :  V 

Bolling  Air  Force  Base  'i  '  ' 

Bldg.  410  ^ 

Washingotn,  DC  20332-6448 

10.  SPONSORING  ,' MONITORING 

AGENCY  REPORT  NUMBER 

AFOSR  90-0039 

11.  SUPPLEMENTARY  NOTES  ftaii  W 

i  1  fe  n  i  ' 

k  ;  A 

-■■  A.rr  T..  1.,  .  .,  n  ,..; 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  'VA  l\P^  A  , 

APPROVED  FOR  PUBLIC  RELEASE:  /.  , 

DISTRIBUTION  UNLIMITED 

^■5 

12b.  DISTRIBUTION  CODE 

UL 

13.  ABSTRACT  (Maximum  200  words) 

This  research  program  involves  the  basic  understanding  and  applications  of 
nonlinear  wave  propagation.  A  number  of  significant  results  have  been  obtained. 
Contributions  include  the  following.  The  analysis  and  computation  of  solitons  and 
their  interaction  properties  in  nonlinear  optical  fibers;  solutions  and  general 
behavior  of  a  class  of  physically  significant  multidimensional  and  novel 
nonlinear  wave  equations;  studies  of  discrete  nonlinear  evolution  systems  and 
numerically  induced  chaos  in  a  class  of  integrable  nonlinear  systems.  This  report 
also  includes  a  list  of  recent  publications  and  preprints. 

19950331  11T 


.  wrMVl> 

Nonlinear  Wave  Equations,  Solitons  and 

Analys is 

IS.  NUMBER  OF  PAGES 

7 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

UNCLASSIFIEI 

1 

UNCLASSIFIED 

UNCLASSIFIED 

SAR 

Standard  -orm  :98  ;9ev  2-89) 
ov  '’o  i 

t’98-^02 


MSN  75^0*0l-280-5500 


FINAL  TECHNICAL  REPORT  (per  AN) 
1  DEC  89-30  NOV  93 


Air  Force  Office  of  Scientific  Research 

"Nonlinear  Wave  Propagation" 

USAF  90-0039 

Principal  Investigator:  Mark  J.  Ablowitz  j 

f 

i 
} 

Program  in  Applied  Mathematics  j 

University  of  Colorado  -  Boulder  i 

Boulder,  CO  80309-0526  ! 

December  22,  1993  ^ 

Abstract 

This  research  program  involves  the  basic  understanding  and  applications  of  nonlinear  wave 
propagation.  A  number  of  significant  results  have  been  obtained  Contributions  include  the 
following.  The  analysis  and  computation  of  solitons  and  their  interaction  properties  in  nonlinear 
optical  fibers;  solutions  and  general  behavior  of  a  class  of  physically  significant  multidimensional 
and  novel  nonlinear  wave  equations;  studies  of  discrete  nonlinear  evolution  systems  and 
numerically  induced  chaos  in  a  class  of  integrable  nonlinear  systems.  This  report  also  includes  a 
list  of  recent  publications  and  preprints. 


1 


1 


Summary  of  Research  Activities 

During  the  past  two  years  there  have  been  eight  journal  publications,  three  more  have  been 
accepted  for  publication  in  journals,  seven  book  chapters  have  been  pubhshed  and  a  monograph 
"Solitons,  Nonlinear  Evolution  Equations  and  Inverse  Scattering"  was  published  by  Cambridge 
University  Press.  In  addition  there  are  nine  Program  in  Applied  Mathematics  Reports  (preprints  of 
the  Program)  which  have  been  submitted  for  publication  and  the  PI  has  been  invited  to  give 
fourteen  lectures  at  universities  and  international  conferences  throughout  the  world.  The  past  two 
years  has  been  a  very  active  and  productive  period  in  this  research  program.  A  brief  overview  and 
description  of  our  research  investigations  are  given  below.  Full  details  are  included  in  the 
publications. 

Soliton  Propagation  in  Nonlinear  Optics 

In  the  application  of  optical  data  communication,  linear  transmission  systems  woric  well 
over  short  distances  or  with  relatively  low  data  rates.  However,  chromatic  dispersion,  dissipation 
and  nonlinear  effects  seriously  limit  the  distances  over  which  the  linear  transmission  is  viable.  In 
practice,  numerous  electronic  repeaters  to  regenerate  linear  signals  are  required  in  order  to  make 
long  distance  data  communication  effective.  Using  nearly  linear  waves  has  a  number  of 
disadvantages  and  researchers  have  studied  the  possibility  of  developing  practical  transmission 
systems  which  use  the  underlying  nonlinearity  of  the  optical  fiber  as  a  critical  feature.  The  net 
result  is  that  solitons  have  proven  to  be  an  important  technological  mechanism  for  data 
communication.  Much  of  the  research  has  centered  on  sending  widely  separated  soUtons  with 
many  equally  spaced  amplifiers  to  counteract  dissipation.  The  amplifiers  are  considerably  easier  to 
produce  than  electronic  repeaters  and  being  optical  in  nature  they  have  the  potential  of  allowing 
nonlinear  systems  to  transmit  information  at  much  larger  data  rates  than  those  of  linear  systems. 

In  many  of  the  important  applications,  the  dominant  underlying  nonlinear  equation  is  one  of 
the  integrable  nonlinear  systems  which  the  PI  and  colleagues  have  intensively  studied  with  the 
support  of  the  Air  Force.  In  the  mathematical  formulation,  solitons  are  associated  with  eigenvalues 
of  certain  linear  systems  and  without  damping  or  amplification,  the  eigenvalues  are  preserved  in 
time.  In  fact  even  with  perturbations,  it  is  known  that  these  eigenvalues  only  evolve  slowly  in 
time.  One  of  the  realistic  models  governing  singly  polarized  optical  fibers  is  the  nonlinear 
Schrodinger  equation,(NLS).  With  the  NLS  equation,  the  real  part  of  the  eigenvalue  is  related  to 
the  carrier  frequency  of  the  optical  wave  and  the  unaginary  part  of  the  eigenvalue  is  associated  with 
the  amplitude.  In  our  research  we  are  interested  in  wavelength  division  multiplexing  (WDM), 
which  is  the  study  of  sending  many  solitons  in  a  given  window  of  time.  The  goal  is  to  transmit  a 
group  of  solitons  which  can  be  measured  easily  at  the  receiving  station.  In  other  words,  we  are 
really  interested  in  transmitting  and  receiving  a  collection  of  eigenvalues  over  long  distances. 
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In  this  project  we  have  collaborated  with  a  group  led  by  Professor  Jon  Sauer  of  the 
Optoelectronics  Computing  Center  at  the  University  of  Colorado.  Mathematically  and  physically 
the  problems  that  must  be  addressed  are  considerably  more  difficult  when  multisolitons  are 
present.  One  difficulty  is  the  need  to  distinguish  separate  solitons  (eigenvalues)  when  they  overlap 
in  physical  space.  In  our  research,  we  have  shown  both  analytically  and  numerically  that  for  two 
solitons  the  power  spectra  (amplitudes  of  the  Fourier  transforms)  always  remain  well  separated  so 
long  as  the  carrier  frequencies  (the  imaginary  parts  of  the  eigenvalues)  of  the  solitons  are  well 
separated.  These  calculations  pave  the  way  to  understand  N  soliton  WDM  and  provides  a  basis  for 
how  a  suitable  perturbation  theory  with  damping  and  frequency  effects  could  be  developed. 

Multidimensional  and  Novel  Nonlinear  Wave  Equations 

The  well  known  integrable  nonlinear  wave  equations  of  physical  significance  include,  the 
Korteweg-de Vries  (KdV)  equation  and  it's  2+1  dimensional  analogue:  the  Kadomtsev-Petviashvili 
equation;  the  nonlinear  Schrodinger  equation  and  its  2+1  analogue:  the  Davey-Stewartson  equation, 
and  the  1+1  and  2+1  dimensional  N  wave  interaction  equations.  These  equations  have  been 
extensively  studied  by  us  in  a  variety  of  contexts.  We  are  now  studying  important  discrete 
nonlinear  wave  systems  such  as  the  2+1  Toda  equation  which  reduces  to  the  well  known  Toda 
lattice  equation  in  1+1  dimension.  The  2+1  Toda  equation  admits  a  variety  of  interesting 
phenomena  depending  on  the  choice  of  signs  in  the  governing  equations.  In  the  choice  of  signs 
where  the  initial  value  problem  is  well  posed,  there  are  two  cases,  one  of  which  admits  linearly 
unstable  waves.  The  inverse  scattering  analysis  in  the  unstable  case  results  in  a  coupled  nonlocal 
Riemann-Hilbert  and  DEAR  system.  For  the  other  choice  of  signs,  the  linearized  Toda  equation  is 
ill-posed.  Nevertheless,  interesting  boundary  value  problems  can  be  formulated  and  solved.  In 
one  of  the  subcases  the  boundary  value  problem  requires  that  a  radiation  condition  be  imposed  in 
order  to  obtain  a  unique  solution.  This  radiation  condition  is  an  analogue  of  the  well  known 
Sommerfeld  radiation  condition  required  in  linear  problems,  such  as  the  Helmholz  equation. 

We  have  been  led  to  study  a  nonlinear  differential-delay  Toda  equation,  which  is  a 
reduction  of  the  2+1  Toda  equation.  Differential-delay  problems  arise  frequently  in  physical 
problems.  This  is  the  first  time  we  have  seen  such  an  equation  in  the  field  of  soliton  theory.  In  the 
study  of  the  Toda  differential-delay  equation,  the  analysis  of  the  associated  inverse  scattering 
problem  is  quite  complicated.  The  inverse  scattering  problem  is  novel.  It  requires  one  to  solve  an 
infinitely  coupled  system  of  Riemann-Hilbert  problems  on  adjacent  strips.  In  previous  applications 
the  inverse  problem  involved  solving  a  Riemann-Hilbert  problem  in  half-planes  or  in  sectors.  The 
soliton  solutions  also  arise  in  a  novel  fashion  and  they  are  described  as  bound  states  on  a  suitably 
branched  Riemann  surface. 
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Wave  collapse  has  been  an  area  of  active  study  by  us  as  well.  It  is  well  known  that  the 
multidimensional  (e.g.  2+1)  cubic  nonlinear  Schrddinger  equations  (arising  in  many  areas  of 
physics;  e.g.  nonlinear  optics)  has  solutions  which  blow  up  in  finite  time.  It  is  less  well  known 
that  one  dimensional  nonlinear  (quintic  and  higher  nonlinearity)  analogues  of  the  Korteweg- 
de Vries  equation  also  admit  solutions  which  tend  to  infinity  in  finite  time.  We  have  recently  shown 
that  nonlinear  versions  (cubic  and  higher  nonlinearity)  of  the  Kadomtsev-Petviashvili  (KP) 
equations  possess  similar  blow-up  phenomena.  We  have  studied  this  problem  both  analytically 
and  numerically.  Indeed,  we  have  found  that  in  the  numerical  calculations  one  must  impose 
appropriate  constraints  which  are  inherent  in  KP  systems  and  which  we  had  studied  earlier  in  the 
context  of  the  integrable  KP  equation  (the  quadratically  nonlinear  system). 

Numerical  Chaos— Truncation  and  Roundoff 

In  our  research  we  have  been  studying  modoate  to  long  time  numerical  expentnents 
involving  integrable  systems,  most  notably  the  nonlinear  Schrodinger,  sine-Gordon  and  modified 
Korteweg-de  Vries  equations  with  periodic  boundary  conditions.  In  our  earlier  work,  supported 
by  this  grant,  we  found  that  in  certain  parameter  regimes  associated  with  given  initial  data, 
spurious  joint  spatial-temporal  numerical  chaos  was  generated  by  uuncation  effects  due  to 
nonintegrable  discretizations.  The  phenomena  can  be  traced  to  the  proximity  of  the  initial  data  to 
underlying  homoclinic  manifolds  inherent  in  the  governing  problem.  We  also  developed  a 
numerical  algorithm  which  was  superior,  and  which  is  based  upon  soliton  theory.  The  latter 
scheme  is  referred  to  as  the  integrable  discretization.  For  the  nonintegrable  discretization,  the  chaos 
was  observed  to  disappear  as  the  mesh  was  refined.  In  our  recent  work,  just  published  in  Physical 
Review  Letters,  we  show  that  there  are  parameter  regimes  in  which  roundoff  effects  induce 
spurious  numerical  chaos.  These  roundoff  induced  instabilities  are  associated  even  with  the  best 
numerical  schemes  such  as  the  Fourier  split-step  algorithm  (used  by  many  researchers)  and  the 
integrable  discretization  mentioned  above.  We  have  found  the  analytical  reason  for  this 
phenomenon.  It  turns  out  to  be  due  to  very  small  effects  in  the  eigenvalues  of  the  associated  linear 
scattering  problem  which  governs  the  solution  of  the  original  problem  via  the  Inverse  Scattering 
Transform.  A  subclass  of  eigenvalues  are  indicative  of  homoclinic  structures  in  the  original  PDF, 
and  under  perturbations,  even  on  the  order  of  roundoff,  tiny  changes  in  these  eigenvalues  translate 
into  large  changes  in  the  solution  of  the  PDF.  We  have  carried  out  a  relevant  perturbation  analysis 
of  the  scattering  problem  and  have  shown  that  this  agrees  with  numerical  computations.  Further 
investigations  will  be  published  in  the  near  future.  We  believe  that  the  numerical  instabUities  and 
chaos  mentioned  above  are,  in  fact  typical  of  what  can  occur  when  the  governing  problem  has 
underlying  homoclinic  manifolds  and  will  be  observed  in  a  wide  class  of  Hamiltonian  wave 
systems,  and  not  necessarily  ones  which  are  integrable. 
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1 .  University  of  Kansais,  Mathematics  Department,  “Chaos  in  the  Nonlinear  Schrodmger  Equation: 

Real  or  Fake,”  Jan.  30,  1992;  “On  Reductions  of  the  Self-Dual  Yang-Mills  Equations,”  Jan. 
31,  1992. 

2.  Carlton  University,  Department  of  Mathematics,  Ottowa,  Canada,  “Sohtons,  Computation  and 

Chaos,”  April  11,  1992. 

3.  Kansas  State  University,  Midwest  Geometry  Conference,  “Reductions  of  the  Self-Dual  Yang- 

Mills  Equations  and  New  Integrable  Systems,”  May  1-3,  1992. 

4.  SIAM  Forum  on  AppUed  Mathematics,  Moderator  of  Panel:  “Training  in  Applied  Mathemat¬ 

ics/Preparation  for  Academic  Careers,”  May  15-17,  1992. 

5.  Canadian  Applied  Math.  Society,  University  of  Alberta,  Edmonton,  Canada,  Conference  on 

Wave  Phenomena,  “Nonlinear  Waves,  Integrable  Systems  and  Modular  Forms,”  June  15-18, 
1992. 

6.  NATO  Advanced  Research  Workshop,  Exeter  University,  Exeter  U.K.;  “Chaos  in  the  Nonhnear 

Schrodinger  Equation-  Can  it  be  Avoided?”,  July  14-18,  1992. 

7.  SUNY  Stony  Brook,  Department  of  Mathematics,  “Nonlinear  Waves,  Integrable  Systems  and 

a  Novel  Class  of  ODE’s” ,  December  3,  1992. 

8.  Challenge  of  Teaching  Science  and  Mathematics,  Graduate  Teacher  Program,  University  of 

Golorado  at  Boulder,  “Goherence  and  Chaos”,  January  11,  1993. 

9.  Colorado  School  of  Mines,  Golden,  Colorado,  Sigma  Xi  lecture,  “Solitons  and  Coherent  Struc¬ 

tures”  ,  March  30,  1993. 

10.  International  Conference  on  Inverse  Scattering,  Bad  Honnef,  Germany,  “Inverse  Scattering  and 

Nonhnear  Wave  Equations,  1  -1-  1,  2  +  1  and  Higher  Dimensions”,  May  17-20,  1993. 

11.  NATO  Lectures,  Bilkent  University,  Department  of  Mathematics  and  Faculty  of  Science, 

Ankarra,  Turkey,  “Coherence  and  Chaos — Illustrative  Case  Studies”;  “Reductions  of  Self 
Dual  Yang  Mills  and  New  Special  Functions  in  Integrable  Systems” ,  May  20-25,  1993. 

12.  Conference:  Nonlinear  Evolution  Equation,  Sohtons  and  the  Inverse  Scattering  Transform, 

Oberwolfach,  Germany,  “Numerical  Chaos  Roundoff  Errors  and  Homochnic  Manifolds”,  July 
11-17,  1993. 

13.  FSU-USA  Conference  on  Chaos,  NAS  Center,  Woods  Hole,  Mass.,  “Roundoff  Errors  and 

Homochnic  Manifolds”,  July  19-23,  1993. 

14.  Nonhnear  Optics  Workshop  and  CoUoquium,  Department  of  Mathematics,  University  of  Ari¬ 

zona,  “Roundoff  Errors  and  Homochnic  MEinifolds”,  Sept.  9-11,  1993. 
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